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ABSTRACT 

A  set  of  absorption  coefficients  for  water  vapor  as 
functions  of  temperature  has  been  empirically  dedecued  from 
existing  quantitative  absorption  and  emission  spectra  between 
1  and  22  The  basic  assumption  was  made  that,  for  optically 
thick  gases,  the  curve  of  growth  is  given  by  a  statistical 
model.  The  absorption  coefficients  were  obtained  for  each 
vibration-rotation  band  by  adjusting  the  band  averaged  line 
half  widths  and  line  spacings  such  that  the  integral  of  the 
absorption  coefficients  was  equal  to  the  known  integrated 
intensity  of  the  particular  oand  at  a  given  temperature.  The 
effect  of  foreign  gas  broadening  and  non-resonant  water  molecule 
broadening  was  included.  The  present  results  are  compared  with 
independent  theoretical  and  experimental  results  and  the  agreement 
is  satisfactory . 
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DfERODUCTION 


For  the  purpose  of  analyzing  and  predicting  rocket  plume  radiation, 
it  is  necessary  to  have  available  a  representation  of  the  spectral  prop¬ 
erties  of  all  the  major  radiating  species  as  functions  of  the  relevant 
parameters.  This  representation  should  be  as  simple  as  possible  yet 
detailed  enough  to  allow  for  determining  the  essential  features  of  the 
spectrum. 

The  major  infrared  molecular  radiators  in  present  rocket  exhausts 
are  carbon  dioxide,  water  vapor,  and  several  diatomic  molecules.  The 
spectral  properties  of  carbon  dioxide  and  many  of  the  diatomic  molecules 
have  been  reported  in  the  literature.  However,  there  are  at  present  no 
theoretical  predictions  of  the  complete  water  vapor  spectra  in  the  infrared 
region  at  high  temperatures.  It  is  the  purpose  of  this  paper  to  derive 

-1 

a  set  of  absorption  coefficients  for  water  vapor  between  50  and  7500  cm 
in  the  temperature  range  from  300°  to  3000 °K  from  available  experimental 
data. 

During  the  last  few  years,  a  large  number  of  spectra  of  high  tempera  - 

1-12 

ture  water  vapor  between  1  and  22  microns  have  been  published.  *  Hiese 
spectra  have  been  measured  by  different  techniques  and  cover  the  tempera¬ 
ture  range  from  300°  to  2700°K.  The  optical  depths  range  between  about 
0.2  to  100  cm*atm,  total  pressures  range  between  50  mmHg  and  10  atm,  but 
the  range  of  parameters  covered  at  any  one  temperature  is  rather  limited, 
especially  at  higher  temperatures .  These  measurements  show  that  the 
spectral  emissivity  depends  on  the  following  independent  parameters: 
absorption  coefficient,  wavelength,  temperature,  partial  pressure,  total 
pressure,  line  broadening  ability  of  foreign  gases,  and  pathlength.  The 
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most  important  parameter  is  the  absorption  coefficient  as  a  function  of 
wavelength  and  temperature .  If  sufficient  thin-gas  spectra  at  all  tempera¬ 
tures  were  available,  a  set  of  absorption  coefficients  could  be  derived 
directly.  However,  this  is  not  the  case.  In  order  to  obtain  absorption 
coefficients  from  measured  spectra  in  which  the  gas  was  not  thin,  the 
functional  relationship  between  emissivity  or  absorptivity  and  optical 

depth  (curve  of  growth)  must  be  known.  A  number  of  theoretical  models 

13  l4  i4 

for  the  curves  of  growth  have  been  developed.  *  Goody  has  treated 
the  absorptivity  of  a  gas  whose  spectral  lines  are  distributed  statistically 
both  in  position  and  intensity.  In  the  same  paper  the  applicability  of 
the  statistical  model  to  water  vapor  was  established.  Also,  Howard,  Burch, 
and  Williams^ have  applied  this  model  to  water  at  room  temperature 
and  1  atm  pressure  and  found  good  agreement  over  a  wide  range  of  optical 
depths . 

The  expression  for  the  spectral  curve  of  growth  contains  the  mean 
spectral  absorption  coefficient,  the  optical  depth,  and  a  rotational  fine 
structure  term  which  is  dependent  on  a  collision  parar*eter  and  the  mean 
line  spacing. 

We  propose  to  use  these  curves  of  growth  to  determine  the  mean  spectral 
absorption  coefficients  from  non-thin  water  spectra  by  adjusting  the  fine 
structure  term  in  such  a  way  that  the  integral  of  the  absorption  coefficients 
over  a  given  vibration-rotation  band  results  in  the  known  value  of  the 
band  intensity. 

DETERMINATION  OF  ABSORPTION  COEFFICIENTS  FROM  NON-THIN  GAS  SPECTRA 

In  the  statistical  model  the  spectral  emissivity  is  given  by 

«<■)  =  d) 
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where  W(uj)  is  the  average  equivalent  width  and  d(tu)  is  the  average  line 
14 

spacing.  Goody  has  shown  that  the  curves  of  growth  are  quite  insensitive 

to  the  choice  of  the  line  intensity  distribution  functions.  When  the 

smoothed  absorption  coefficient  and  mean  line  width  to  line  spacing  ratio 

are  chosen  so  that  the  curves  of  growth  agree  in  the  optically  thin  and 

W 

in  the  square  root  limits,  the  values  of  ^  corresponding  to  a  delta  func¬ 
tion,  an  exponential  or  a  l/S  distribution  differ  by  at  most  25$. 

In  the  present  reduction  of  the  experimental  data,  we  have  used  the 
intermediate  form  corresponding  to  an  exponential  distribution  of  line 
intensities : 


where  k(tu)  is  the  mean  absorption  coefficient,  u  the  optical  depth  (at 
STP),  and  a(tu)  the  fine  structure  parameter.  The  fine  structure  term 
a  Ob)  is  proportional  to  the  local  mean  value  of  the  ratio  of  the  collision 
half -width  to  the  line  spacing,  *y/d.  It  is  to  be  noted  that  the  value 
of  d,  deduced  from  experimental  curves  of  growth,  is  only  meaningful  when- 
interpreted  in  terms  of  a  particular  intensity  distribution  function. 

In  the  present  analysis,  we  have  chosen  to  normalize  the  value  of  d 
so  that  it  is  appropriate  for  a  delta  function  distribution  (i.e.,  a  = 

The  line  width  is  expected  to  be  proportional  to  the  number  of 
collisions  experienced  by  the  water  molecule  per  unit  time.  Since 
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collisions  with  other  molecules  of  water  vapor  are  more  effective  than 

17 

those  with  molecules  of  a  different  species,  it  has  been  found  that 

good  agreement  with  the  experimental  data  at  room  temperature  has  been 

obtained  by  assuming  that  the  line  width  is  proportional  to  (Vg^QPg^o  + 

YxPx)>  w*iere  Y  =  line  width  and  the  subscript  x  designates  the  foreign 

gas.  Since  we  want  to  extend  this  treatment  to  higher  temperatures,  the 

effect  of  interactions  between  resonating  and  non-resonating  water  molecules 

lB,19 

and  their  temperature  dependency  should  be  taken  into  account.  Urns, 

we  obtain  a  three  term  expression, 

a(m)  d(x)  =  ^0(^0^  T~  +  ^0^  /  f)  +  Px\^  J  T~  * 

where  y'  is  the  line  width  due  to  non-resonating  dipole  interactions  of 
*2° 

HgO  at  Tq.  The  following  assumptions  are  now  made:  The  line  widths  y °(x) 
due  to  interactions  with  resonating  and  non-resonating  water  molecules  and 
with  foreign  gas  species,  which  are  functions  of  frequency,  can  be  approxi¬ 
mated  by  a  band  averaged  value,  so  that  the  frequency  dependence  of  the 
broadening  term  a(x)  appears  only  through  the  mean  line  spacing  d(u>). 

Thus,  Eq.  (3)  becomes 


a(x)  = 


VH_0  Fr  ,  F  N 

dW  pt  •/  T  [c  (  /  T  *  ”')*  (1'c)ax] 


(4) 


where  a'  =  /0  ' ° 


and  ax  =  ^0  and  C  ls 


the  mole  fraction  of  water. 


o  -  17-20 

Values  of  y„  _  and  a  are  given  in  the  literature  and  the  following 

*2°  X  —  ^1  _i 

values  will  be  used  in  the  present  study:  yu  n  =  0.5  cm  atm  ;  =  0.18; 


H,0 
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aQ  =  0.09*  Since  no  value  for  o'  is  given,  we  will  assume  initially  that 
o'  is  approximately  0.1.  As  a  next  step  we  assume  that  a  band  averaged 
mean  line  spacing  d  can  be  used  so  that  a(tu)  becomes  independent  of  the 
frequency  and  we  can  write  Sq.  (2)  as 


The  frequency- independent  broadening  term  a'  can  then  be  evaluated  from 

w 

existing  spectra,  where  — (tt) )  =  -  0n  (l-e[tt)])>  by  using  as  the  second  of 
the  two  equations  required  to  define  k(aj)  and  a 

J  k(o>)du)  =  aO0  (T ) ,  (6) 

band 

where  a°  is  the  known  integrated  intensity  at  STP  of  the  band  in  question 
and  0(T)  is  a  known  function  of  the  temperature.  It  can  be  shown  that, 
in  the  temperature  range  of  interest  here  (up  to  ~  3000°K),  the  function 
0(T)  =  1  for  fundamental  bands  in  general.  However,  as  pointed  out  recently 
by  W.  S.  Benedict,  the  6.3-m>  fundamental  band  of  HgO  may  not  obey  this  rule 
because  of  strong  rotation-vibration  interaction.  Room  temperature  measure¬ 
ments  of  the  6.3-h  band  indicate  a  value  of  250  (cm  ^  per  cm  at  STP),  while 
measurements  at  temperatures  of  2700 °K  indicate  a  value  of  3&0  (cm  ^  per  cm 
at  STP).  Therefore,  an  empirical  function  0(T)  =  250  +  0.05T  has  been  used 
to  describe  the  temperature  dependence  of  the  6.3-p.  band  between  300 °K  and 
3000 °K.  For  combination  and  overtone  bands,  0(T)  was  given  in  Ref.  21  for 
the  1.87-  and  1.38-p,  bands.  The  values  of  a°  in  (cm~^  per  cm  at  STP)  for 
the  four  major  vibration-rotation  bands  used  in  the  present  analysis  are 
shown  in  Table  I. 
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Table  I:  Integrated  intensities  used  in  present  analyses. 
Bands  6.3  p.  2.7  u  1.87  u  1.38  U 

250  ±  20$  230  ±  155&  26  ±  15i>  21  ±  15* 


The  pure  rotational  band  could  not  be  evaluated  in  this  fashion  since 
no  measurements  for  the  total  band  exist.  However,  mean  absorption  coef¬ 
ficients  were  calculated  in  Ref.  11;  these  have  been  included  in  the 
present  tabulation. 

Approximately  80  published  emission  and  absorption  spectra  were 
used  in  the  present  analysis.  The  spectral  emissivity  at  intervals  of 
25  cm-1,  together  with  the  optical,  depth,  temperature  and  total  pressure 
served  as  inputs  in  a  computer  program.  When  the  experimentally  determined 
spectra  were  given  with  a  higher  resolution  than  25  cm  \  the  nr  »  re 
smoothed  by  hand.  The  absorption  coefficient  was  calculated  from 


by  choosing  a  such  that  Eq.  (6)  was  satisfied  (see  Fig.  l).  If 

-  I  to  1  y  *  du)  «  a  ,  ( 

u  J  l-e(tu)  * 

the  spectrum  was  from  a  thin  gas  sample  and  the  absorption  coefficients 
could  be  determined  directly. 

Two  sources  of  error  in  the  determination  of  a(T)  and  k(u,T)  are 
recognized : 

l)  Error  in  the  measurements  of  the  state  parameters  of  the  gas 
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or  of  the  spectral  emissivity  or  absorptivity,  and 
2)  Uncertainty  in  the  integrated  intensities  of  the  vibration- 
rotation  bands. 

An  analysis  was  performed  in  which  reasonable  error  limits  for  the  state 
parameters  and  the  integrated  intensities  were  selected.  Hie  state  param¬ 
eters  T,  u,  c,  p^,  and  the  emissivity  were  varied  by  ±  5$,  while  or  was 
varied  by  ±  20$.  It  is  believed  that  these  variations  constitute  upper 
limits,  since  many  experimenters  quote  smaller  error  limits.  It  is 
found  that  the  uncertainty  of  ±  3$  in  the  measured  emissivity  has  the 
greatest  influence  on  a(T)  and  k(<u,T)  (Aa  s=  ±  100$,  Ak  ±  20$).  An 
uncertainty  of  ±  20$  in  a  gives  Aa  se  ±  20$  and  Ak  a  i  15$.  All  other 
parameters  influence  a  and  k  by  ±  5$  or  less . 

RESULTS 

The  average  absorption  coefficients  obtained  by  the  method  described 
in  the  previous  section  are  tabulated  in  Table  II,  and  are  plotted  in  Fig.  2 
for  the  two  extreme  temperatures,  300°K  and  3000 °K.  In  general,  the  mean 
deviation  of  the  absorption  coefficients  taken  from  the  individual  spectra 
is  within  ±  20$.  As  an  example,  in  Fig.  3  k  at  1250  cm-1  is  plotted  versus 
temperature.  The  solid  line  represents  the  mean  of  the  individual  points. 
From  this  curve,  the  absorption  coefficients  at  3008,  600°,  1000°,  1500°, 
2000°,  2500°,  and  3000°K  are  selected  and  entered  in  Table  II.  In  some 
portions  of  the  spectra,  where  the  absorption  coefficient  changes  rapidly, 
a  much  greater  spread  in  the  individual  absorption  coefficients  is 
observed.  In  Fig.  ^  k  at  3750  cm  \  the  Q-branch  region  of  the  2. 7-p. 
band,  is  plotted  versus  temperature.  The  large  spread  here  is  probably 
introduced  by  small  errors  in  the  wavenumber  calibration.  Also,  the 
absorption  coefficients  in  the  spectral  regions  between  the  individual 
bands  ("troughs”)  and  in  the  region  tu  >  7500  cm  ^  are  not  very  certain 
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because  of  lack  of  sufficient  data.  Since,  at  the  present  time,  there 
are  no  data  for  temperatures  below  2000 °K  in  this  region,  no  absorption 
coefficients  are  listed  in  Table  II. 

The  average  line  s pacings  d(T)  obtained  from  the  values  of  a(T), 
using  Eq.  (4),  are  plotted  in  Fig.  5*  Data  from  approximately  80  spectra 
of  the  2.7~  and  6.3-p  band  are  represented.  Each  point  has  about  a  factor 
of  2  uncertainty.  The  overall  spread,  of  points  is  mostly  within  a  factor 
of  2  to  3,  but  can  be  as  large  as  a  factor  of  10  (at  1200 °K  for  instance). 
The  straight  line  is  a  least  square  fit  to  the  points  and  is  given  by 

d(T)  =  exp[-  .00106T  +  1.21]  (9) 

One  point  at  2200 °K  is  obtained  from  the  rotational  band.  The  point  from 
Oppenheim  ^  at  1200°K  is  a  simple  average  taken  from  his  measured  (2ny  /d) 
values  (see  Fig.  6).  The  majority  of  points  in  Fig.  5  are  restricted  to 
temperatures  <  l800°K  because  the  high-temperature  measurements  made  with 
the  rocket  motor  in  this  laboratory  were  in  the  thin  gas  region  where 
Eq.  (8)  applies. 

The  two  combination  bands  at  1.87  ^  and  I.38  p  are  much  weaker  than 

Q 

the  fundamental  bands.  Room  temperature  absorption  spectra  of  the  1.87-p 
band  indicate  that  the  value  of  d  is  practically  the  same  as  that  obtained 

from  the  two  fundamental  bands .  Values  of  d  derived  from  the  data  presented 

3  — 

by  Nelson  are  generally  higher  than  values  of  d  obtained  from  his  data 

for  the  two  fundamental  bands,  but  the  experimental  uncertainty  in  his 
measurements  does  not  permit  any  definite  conclusions.  Our  own  data  were 
always  obtained  at  low  optical  depths,  so  that  no  d  could  be  obtained. 

In  the  absence  of  further  information,  it  is  assumed  that  the  d(T)  for 
the  two  combination  bands  is  the  same  as  that  deduced  for  the  fundamental 
bands . 
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The  rotational  band  is  very  strong  and  fine  structure  affects  the 
measurements  even  at  the  highest  temperatures  measured  (~  2700°K).11 
However,  since  the  measurements  were  carried  out  only  up  to  22  n,  no 
total  band  strength  could  be  obtained.  Thus,  the  procedure  used  before 
for  the  vibration-rotation  bands  could  not  be  employed.  Additional  know¬ 
ledge  about  the  integrated  intensity  of  this  band  up  to  22  n  was  needed 
which  presupposed  the  knowledge  of  the  spectral  absorption  coefficient. 

The  theoretical  relationships  of  Ref.  11  were  used  here  in  a  semi -empirical 
fashion  by  adjusting  the  effective  mean  rotational  constant  slightly  to  fit 
the  observed  spectra  in  the  far  wing  of  the  rotational  band.  The  values  of 
k  thus  determined  were  used  to  calculate 

1000  cm"1 

[  k(uj,T)  du>  =  c/(T)  (10) 

450  cm-1 

which  then  replaced  Eq.  (j). 

COMPARISON  WITH  OTHER  THEORETICAL  AND  EXFERTMEJfTAL  RESULTS 

A  check  of  the  absorption  coefficients  derived  here  against  indepen- 

24  25 

dent  work  is  possible  for  some  selected  temperatures.  Benedict  *  has 
obtained  absorption  coefficients  for  individual  lines  in  the  6.3-n  band 
at  «s  300 °K  and  as  600°K  and  in  the  2.7-^  band  for  as  300 °K.  We  have  averaged 
Benedict’s  values  and  compared  them  with  the  absorption  coefficients  taken 
from  this  work  (see  Figs.  7,  8,  and  9)-  Oppenheim’s  detailed  measure¬ 
ment  of  the  absorption  coefficients  in  the  2.7-H  band  at  1200 °K  is 

7 

compared  in  Fig.  10.  Goldstein’s  determination  of  the  absorption  coef¬ 
ficients  of  the  I.38-,  I.87-,  2.7-'.,  and  6.3-n  bands  at  as  500  °K  and  1200 °K 
are  compared  with  the  present  values  in  Figs.  11  through  l6.  One  set  of 
experimental  data  has  recently  become  available.  Simmons,  Arnold,  and 
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Smith  have  measured  the  2.7-n  bands  between  700°  and  1200 °K  at  various 
optical  paths.  Since  their  data  were  not  included  in  our  derivation  of 
the  absorption  coefficients  and  band  averaged  line  spacing,  they  provide 
us  with  an  independent  check  for  calculating  thick  gas  spectra  of  the 
2.7-u  bands  at  around  1000°K.  The  result  of  this  comparison  for  some  of 
the  spectra  is  given  in  Figs.  17-22.  Similar  agreement  is  found  with  the 
other  spectra  from  Ref.  26. 

In  a  consistency  check,  some  of  the  measured  spectra,  which  were 
used  to  obtain  a  band-averaged  line  spacing  and  average  absorption 
coefficients,  are  recalculated  using  the  d  from  Fig.  5  (dashed  line) 
and  the  k's  from  Table  II.  In  Fig.  23  we  compare  Nelson's  measure¬ 
ments  of  the  6.3-u  band  at  1111 °K  with  the  present  representation  of 
k (uj )  and  T.  In  Figs.  24,  25,  and  26,  the  spectra  obtained  in  this 

laboratory  are  compared  at  2200°K.  In  Fig.  27,  one  spectrum  obtained 

2 

by  Burch  and  Gryvnak  of  the  2.7-U  band  at  1200 °K  is  compared  with  the 
present  results.  In  general  the  agreement  is  adequate  for  the  present 
representation  of  H^O  at  high  temperatures . 

CONCLUDING  REMARKS 

The  absorption  coefficients  averaged  over  25  cm  1  intervals  given 
in  this  study  for  temperatures  greater  than  300°K  represent  the  'best" 
values  available  at  the  present  time.  (The  values  at  300 °K  were  included 
here  only  for  completeness.  More  detailed  tabulations  are  available  in 
the  literature,  for  example,  References  2k  and  27).  The  absorption 
coefficients  derived  here  have  been  based  on  nearly  all  the  available 
published  spectra,  using  a  statistical  model  for  the  curve  of  growth 
and  a  frequency-independent  fine  structure  term  to  extrapolate  to  zero 
pathlength.  In  the  slowly  varying  portions  of  the  spectra,  the  uncer- 
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tainties  are  believed  to  be  within  ±  20^,  while  in  portions  of  the  spectra 
where  steep  changes  occur,  the  uncertainties  may  be  higher.  Also  the 
values  in  the  troughs  between  the  vibration- rotation  bands  are  not  very 
certain  because  of  the  lack  of  sufficient  data. 

The  different  frequency-independent  fine  structure  terms  obtained 

from  the  analysis  of  the  optically  thick  spectra  were  used  to  calculate 

band-averaged  line  s pacings  as  functions  of  temperature.  These  individual 

points  are  uncertain  by  a  factor  of  2  and  although  the  spread  can  be 

as  high  as  a  factor  of  10  (at  1200°K  for  instance),  the  mean  values 

(dashed  curve  in  Fig.  5)  are  believed  to  be  uncertain  by  a  factor  of  2 

only.  Biat  these  mean  values  may  be  used  for  first  approximations  in 

the  calculation  of  moderately  thick  gases  was  shown  in  the  comparison 

26 

with  some  of  the  spectra  by  Simmons,  et.al.,  in  Figs.  17-22.  More 
accurate  calculations  for  very  thick  gases  will  require  a  detailed 
evaluation  of  the  dependence  of  the  fine  s tincture  parameters  on  both 
temperature  and  wavelength. 
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ntegrated  intensity  a  of  the  2.7-d  band  as  function  of  a  for  various  temperatures 


001 


Figure  2.  Plot  of  absorption  coefficients  at  Z00°K  and  3000*K 


T  (°K) 


Figure  3-  Compilation  of  average  absorption  coefficients  versus  temper¬ 
ature  at  1250  cm”^  as  obtained  from  spectra  taken  by  several 
investigators.  Solid  line  is  the  average  value 
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Oppenheim  of  the  2.7-U-  band  at  1200°K. 
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Figure  7* 
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Comparison  of  absorption  coefficients,  obtained  by  Benedict 
at  300 *K,  averaged  over  25  cm  \ 
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Figure  6.  Comparison  of  absorption  coefficients,  obtained  by  Benedict 
at  600%  averaged  over  25  cm  \ 
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Figure  10.  Comparison  of  absorption  coefficients,  obtained  by  Oppenhe 
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Figure  11.  Comparison  of  absorption  coefficients .  obtained  by  Goldstein 
of  the  6.3-ix  band  at  i4-73°K  (solid  line).  Present  values  are 
given  as  points . 
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Figure  13.  Comparison  of  absorption  coefficients,  obtained  by  Goldstein 
of  the  1.9-p  band  at  U73°K  (solid  line).  Present  values  are 
given  as  points . 


26 


6800  6900  7000  7100  7200  7300  7400  7500  7600 

WAVENUMBER  (CM-1) 


WAVENUMBER  (CM' 


30Nvidaosav 


32 


Figure  19.  Comparison  of  measured  and  calculated  emissivity  of  the  2.7-U  band.  Solid  lines  are  t.he 
measured  spectra  taken  from  Ref.  26,  points  are  the  present  values. 
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the  measured  spectra  taken  from  Ref.  26,  points  are  the  present  values. 


Figure  21.  Comparison  of  measured  and  calculated  emissivity  of  the  2.7-u  band.  S 
the  measured  spectra  taken  from  Ref.  26,  points  are  the  present  values 
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ABSORPTION  COEFFICIENTS  OF  H20  (CONT.O) 
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ABSORPTION  COEFFICIENTS  OF  H20  (CONT.O) 
K  (PER  CM  AT  STP ) 
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ABSORPTION  COEFFICIENTS  OF  H20  (CONT.D) 

K  (PER  CM  AT  STP ) 
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